      function cmstar(xmu)

      implicit real*8(a-h,o-z)
      cmstar=6.8741+2.6933*xmu+.3601*xmu*xmu-.0154*xmu*xmu*xmu
      return
      end